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Abstract 

The nonlocal properties for a kind of generic N- dimensional bipartite 
quantum systems are investigated. A complete set of invariants under 
local unitary transformations is presented. It is shown that two generic 
density matrices are locally equivalent if and only if all these invariants 
have equal values in these density matrices. 
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As a fundamental phenomenon in quantum mechanics, nonlocality has been given a lot 
of attention in foundational considerations, in the discussion of Bell type inequalities [1] 
and hidden variable models, see e.g. [2]. Nonlocality turns out to be also very important 
in quantum information processing such as quantum computation [3], quantum teleporta- 
tion [4, 5, 6, 7], dense coding [8] and quantum cryptographic schemes [9, 10, 11]. Nonlocal 
correlations in quantum systems imply a kind of entanglement among the quantum sub- 
systems. The nonlocal properties as well as the entanglement of two parts of a quantum 
system remain invariant under local transformations of these parts. 

The method developed in [12, 13], in principle, allows one to compute all the invariants 
of local unitary transformations, though it is not easy to perform it operationally. In [14] , 
the invariants for general two-qubit systems are studied and a complete set of 18 polynomial 
invariants is presented. It is proven that two qubit mixed states are locally equivalent if 
and only if all these 18 invariants have equal values in these states. Therefore any nonlocal 
characteristics of entanglement is a function of these invariants. In [15] three qubits states 
are also discussed in detail from a similar point of view. 

In the present paper we discuss the locally invariant properties of arbitrary dimensional 
bipartite quantum systems. We present a complete set of invariants and show that two 
generic density matrices with full rank are locally equivalent if and only if all these invariants 
have equal values in these density matrices. 

We first consider the case of pure states. Let H be an iV-dimensional complex Hilbert 
space, with \i), i = 1, N, as an orthonormal basis. A general pure state on H <g> H is of 
the form 

N 

I*) = E %K> ® li>. a a e € (!) 

N 

with the normalization ^ a^a*,- = 1 (* denoting complex conjugation). 

A quantity is called an invariant associated with the state |^) if it is invariant under all 
local unitary transformations, i.e. all maps of the form U ®U from H <g> H to itself, where 
U is a unitary transformation on the Hilbert space H. Let A denote the matrix given by 
(A)ij = ciij. The following quantities are known to be invariants associated with the state 
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given by (1), see [16, 17, 18, 19]: 

I a = Tr(AAY, a = l,...,N; (2) 

(with the adjoint of the matrix A). 

In terms of the Schmidt decomposition, a given can always be written in the following 
form, using two suitable orthonormal basis i = 1, ...,iV: 

AT 

l*> = EV^I < > / ' 



where ^ Aj = 1, Aj > 0. The Aj, i — 1, A/", are the eigenvalues of the matrix AA^ . As 

i=i 

AA^ is self-adjoint, there always exists a unitary matrix V, VV^ = VW = 1, such that 
VAAW = diag{A u ...,A N }. The invariants (2) can then be written in the form: 

N 

/ Q = £A?, a=l,...,N. 
i=i 

As the eigenvalues of the matrix are given by the invariants under local unitary 
transformations, two pure states (on H <g> H) are equivalent under locally unitary trans- 
formations if and only if they have the same values of the invariants I a , a = 1, ...,N [20]. 
Moreover two Hermitian m x m matrices A and B are unitary equivalent {i.e., there exists 
a unitary matrix u on an m-dimensional complex vector space satisfying A = uBv)) if and 
only if 

Tr(A a ) =Tr(B a ), for a = 1, m. (3) 

We consider now mixed states on H ® H. Let p be a density matrix defined on H ® H 
with rank(p) = n < N 2 . p can be decomposed according to its eigenvalues and eigenvectors: 

n 

p = X^NK^I' 
i=i 

where Aj and |z/j), z = 1, ...,n, are the nonzero eigenvalues and eigenvectors respectively of 
the density matrix p. |z/j) has the form 

N N 

Wi) = E 4i|fc>® 10. a Le<L7, ^44 = 1, i = l,...,n. 
fc,/=i fe,/=i 

Let denote the matrix given by (Ai) k i = a l kl . We introduce {pi}, {#;}, 

Pi = Tr 2 \vi){vi\ = AiA\, 6 i = (Tr 1 \v t ){v l \)* = A\A^ i, j = 1, ...,n, (4) 



Tri and Tr 2 stand for the traces over the first and second Hilbert spaces respectively, and 
therefore, pi and ^ can be regarded as reduced density matrices. Let and @(p) be two 
"metric tensor" matrices, with entries given by 

ft(pk = TV(p iPj ), e(p) -=Tr(^), fori,j = l,..,n, (5) 

and 

to(p)ij = ®(p)ij = 0, for iV 2 > i,j > n. 

We call a mixed state p a generic one 5 if the corresponding "metric tensor" matrices Q, 6 
satisfy 

det(Q(p)) ^ 0, and det{Q{p)) ^ 0. (6) 

Obviously, a generic state implies n = N 2 or det(p) ^ 0, namely, a state with full rank. 
Nevertheless a fully ranked density matrix may be not generic in the sense of (6). 

Similarly we also introduce X(p) and Y(p) as 

X(p) ijk = Tr(p i p j p k ), Y(p) ijk = Tr(9 i 6 j 6 k ), i, j, k = 1, n. (7) 

[Theorem]. Two generic density matrices with full rank are equivalent under local unitary 
transformations if and only if there exists a ordering of the corresponding eigenstates such 
that the following invariants have the same values for both density matrices: 

J s (p)=Tr 2 (Tr lP s ), s = l,...,iV 2 ; 

(8) 

Q(p), 6(p), X(p), Y(p). 

Remark 1. It is well-known that the set of eigenvalues and corresponding eigenstates is 

uniquely defined, but not the labelling of them. However, from the proof below, one can 

see that two generic density matrices would have the same set of eigenvalues if they share 

the same values {J s (p)}. One can uniquely choose the label for the eigenstates with the 

different eigenvalues. For the case of degenerate eigenvalues, if two generic density matrices 

p and p' are equivalent under the local unitary transformations, one can always find a kind 

of label for the eigenstates such that they share the same invariants (8), i.e., under this 

label, n(p)ij = to(f/)ij, &(p)ij = &(p')ij, X(p) ijk = X(p') ijk , Y(p) ijk = Y(p') ijk . This is due 
5 These states are all the ones but a set of measure zero: {p \det(fl(p)) — 0, det(Q(p)) = 0}. 
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to that these invariants are the sufficient and necessary conditions for two generic density 
matrices to be equivalent under local unitary transformations, see the proof below. 

[Proof]. We first show that the quantities given in (8) are invariant under local unitary 
transformations. Let u and w be unitary transformations, uv) = v)u = ww^ = w^w = 
1. Under the local unitary transformation u (g) w, we have p — > p' — u <g> w p <g> . 
Correspondingly, we have \vi) — > = u <g> w\fi), or equivalently Ai is mapped to A\ = 
vtAiW, where u l is the transpose of u. Therefore 

p\ = A\A} = u'AAju* = u l Pl u\ 9{ = A'}A\ = w^A\A lW = w^ iW . (9) 

By using (9), it is straightforward to check the following relations: 

n n 

r{ P ) J S ( P >) = tv 2 e AjTnCMM)] = rr 2 E = J s (p), 

1=1 1=1 
nip)^ - n{p% = Tr{p' lP ' 3 ) = Triu'p^u*) = Q(p) lv 

Q(P)H "> ®(P% = TriO'^) = Tr(w%e jW ) = 9(p) y 

X(p) ijk -> X(p') ijk = Trirtp'jPk) = Triii* piPjp k u*) = X(p) ijk , 

Y(p) ijk - Y(p% k = Tr{Q\d%) = Tr(w%6je k w) = Y(p) tjk , 

where k = 1, ...,n. Hence the quantities in (8) are invariants of local unitary transfor- 
mations. If two density matrices are equivalent under local unitary transformations, then 
their corresponding invariants in (8) have the same values. 

Now suppose conversely that the states p = Z)"=i -M^K^I an d p' = Yh=i KWDi^ll & ve 
the same values to the invariants in (8). We are going to prove that p and p' are equivalent 
under local unitary transformations. 

a) As 

J S {P) = Tr 2 (±\tT ri (\u z )(u t \)) = Tr 2 (±XtAAl) = ±\°, 

i=l i=l i=l 

from J s (p') = J s {p) we have 

n n 

£A? = EAJ> V S = 1,...,A 2 . 

i=l i=\ 

From (3), we have that p' and p have the same nonzero eigenvalues, i.e., A^ = Aj, i — 1, n. 

b) From (5), the generic condition det(Q(p)) ^ implies that {pi}, % — 1, ...,n(= A 7 " 2 ), 
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span the space of N x N matrices and 

n 

PiPi = T, c vPk> c l e c ( 10 ) 
k=i 

Taking trace of (10) and using the condition Trpi = 1 one gets 



k=l 



From (11) and (7) we obtain 



Therefore 



i=i 



cj,- = x;^*n ,fc . (i 2 ) 
fc=i 

where the matrices fi y is the corresponding inverses of the matrices Qij (which exist due 
to the assumption (6)). (12) implies that the coefficients C\j are given by {f2y, X ijk }. 
From (5), the generic condition (6) implies that {p^ forms an irreducible iV-dimensional 
representation of the algebra gl(N, C) with the generators {e i: i — 1, A^ 2 } satisfying 

N 2 

M^ffcx' ( 13 ) 

k=l 

where = C^ — C^. More explicitly, it (e,) = pi, i — 1, A^ 2 , where 7r is the representation 
of gl(N, C). 

The generic condition det(Q(p J )) ^ implies that {p-}, z = 1, A^ 2 , also span the space 
of N x N matrices, 

n 

pW = E<%V*, Cy^C. ( 14 ) 
fe=i 

If ft(p') = Q(p) and X(p') = X(p), we have Cg = C^-. Therefore, { Pi } and {p^} (if one 
chooses 7r'(ej) = p-) are two irreducible A^-dimensional representation of g>/(A^, C) (13). It is 
well-known that all the Casimir operators of the algebra gl(N, C) can be expressed in terms 
of homogeneous polynomials of e^'s (for example, the first Casimir operator C 2 can be written 
as a quadratic polynomial of e^'s). Moreover Casimir operators are algebraically independent 
and give rise to a complete set of generators for the center of the universal enveloping algebra 
of gl(N, €). They take scalar values on an irreducible representation of gl(N, C) (from 
Schur's Lemma), and become the characters of the irreducible representations [21]. Due to 
the fact that the trace of every polynomial of {p^} and {p^} can be represented in terms 
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of {flij(p), Xij k (p)} and {fiy(p'), Y^^p')} respectively (see below remark 2), we conclude 
that the values of all the Casimir operators given by the two N- dimensional representations 
{pi} and {p'j} are equal, from the condition f2(p) = fl(p') and Xijk(p) = Xij k (p'). Hence, the 
two sets of representations (primed and unprimed) of the algebra gl(N, C) are equivalent, 
i.e., 

p\ = u'piU*, (15) 

for some u eW. 

Similarly, from 0(p) = 0(p') and Y ijk (p) = Y^- fc (p') we can deduce that 

9[ = w^Oiw, for some w G U. (16) 

From the Singular value decomposition of matrices [22], we have |z/) = u <g) w\i/i), i = 
1, A^ 2 , and p' = u ® w p ® w^. Hence p' and p are equivalent under local unitary 
transformations. ■ 

Remark 2. For a degenerate state p, det(£l(p)) = (resp. det(Q(p)) = 0), the above 
invariants (8) are not complete in the sense that two degenerate density matrices can be 
not equivalent under local unitary transformations even if they give the same values to the 
invariants in (8). This is due to the fact that there exist null vectors for the degenerate 
state. For example in the case det(Q(p)) = 0, there exists at least one Hermitian matrix 
B which satisfies Tr(Bpi) = for i = l,...,n. Hence fi(p')ij = Q(p)ij and X(p')ij k = 
X(p)ij k are not enough to get the first equivalence relation (15). In this case some new 
invariants have to be introduced to get a complete set of invariants. From the algebraic 
relations (10) and formula (12), other generalized invariants like Tr((pi) mi (pj) m: > ...(p k ) mk ) 
and Tr((9i) mi (9j) m: > ...(9 k ) mk ), i,j, k — 1, n; m^, rrij, m k G IV, can be represented in 
terms of {flij, Xij k } and {0jj, Yij k } for a generic state with full rank, for example, 

Tr(n- n- . . . n \ — r ai C 012 . . . n am - 2 O 

Hence by doing so we do not get new invariants. 

To summarize, we have discussed here the local invariants for arbitrary dimensional 
bipartite quantum systems and have presented a set of invariants of local unitary transfor- 
mations. The set of invariants is not necessarily independent (they could be represented by 
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each other in some cases) but it is complete in the sense that two generic density matrices 
are equivalent under local unitary transformations if and only if all these invariants have 
equal values for these density matrices. 
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